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ABSTRACT 

Kapteyn’s star is an old M subdwarf believed to be a member of the Galactic halo population of 
stars. A recent study has claimed the existence of two super-Earth planets around the star based 
on radial velocity (RV) observations. The innermost of these candidate planets-Kapteyn b (P = 48 
days)-resides within the circumstellar habitable zone. Given recent progress in understanding the 
impact of stellar activity in detecting planetary signals, we have analyzed the observed HARPS data 
for signatures of stellar activity. We find that while Kapteyn’s star is photometrically very stable, 
a suite of spectral activity indices reveals a large-amplitude rotation signal, and we determine the 
stellar rotation period to be 143 days. The spectral activity tracers are strongly correlated with the 
purported RV signal of “planet b,” and the 48-day period is an integer fraction (1/3) of the stellar 
rotation period. We conclude that Kapteyn b is not a planet in the Habitable Zone, but an artifact 
of stellar activity. 


1. INTRODUCTION 


M dwarfs are very attractive targets in the search for 
potentially habitable terrestrial planets. Their compara¬ 
tively small masses and radii make the de tection of low- 
m ass pla nets in the habitable zone (HZ; [Kopparapu et 


al. 2013) via radial velocity (RV) and transit photome¬ 


try feasible even now. The observational advantages of 
M dwarf planets also makes them am enable to atmo¬ 


spher ic characterization (e.g. GJ 1214b; Kreidberg et al. 


2014), and there is a possibility that JWS'i' could place 


constraints on the atmospheric composition-including 
biosignat ures-of a terrestrial planet in the HZ of a nearby 


M dwarf ( Seager et al.||2013 Cowan et al.||2015 ). 


An insidious problem for RV detection of low mass 
planets around M dwarfs, however, is the presence of RV 
signals created by stellar activity. Starspots, plage, fil¬ 
aments, convective suppression, and activity cycles can 
all create apparent RV shifts by altering the balance be¬ 
tween re d- and blueshifted por tions of the stellar sur¬ 
face (e.g. Dumusque et al.|[2014 and references therein). 
Activity-induced RV signals appear preferentially at the 
stellar rota tion period and its integer fractions 


loisse 


nods around lUU days 1 

Engle & Guinan 

2011 

), periods 

associated with the HZ will be plagued 

oy actiyity sig- 

nals. We have recently shown examples (Robertson et al. 


odicities believed to be low-mass planets in or near the 
HZ were actua lly caused by activ ity. Eurther more, for 
Barnard’s star ( Kiirster et al.|20d3 ) and GJ 581 (Robert¬ 


son et al. 2014), the activity signals which drove RV shifts 


did not create a detectable photometric periodicity, sug¬ 
gesting that certain magnetic phenomena-probably mag¬ 
netic suppression of convection-may only be apparent in 
spectral activity tracers. Magnetic activity signals ap¬ 


pear to be ubiquitous at the RV < 1 m/s level. 

Kapteyn’s star (= GJ 191) is a halo star with high 
proper motion, which may originate from ou tside the 
Galaxy (see review by Kotoneva et al. 2005). At 3.9 
pc from the Sun, it is the nearest halo star. The star 
is a target of the HAR PS exoplanet search program, 
and Bonfils et al. (2013) reported preliminary solutions 
of possible planet d etections. More recently, [Anglada^ 


Escude et al. (2014 hereafter A14) reported the detec¬ 
tion of two super-Earth exoplanets at P = 48.6 (planet 
b) and P = 121.5 (planet c) days, based partially on 
those same HARPS observations. The proposed orbit for 
planet b places it in the HZ. Assuming a stellar age typ¬ 
ical of halo stars, the existence of an HZ planet around 
Kapteyn’s star allows for the fascinating possibility of 
life nearly as old as the Universe itself. Gonsidering the 
ramifications of such a discovery, it is doubly important 
to ensure Kapteyn b truly exists. 

A14 considered the possibility that activity creates sig¬ 
nals in the RVs of Kapteyn’s star. However, two of 
the activity tracers they examined-the Ca H H&K lines 
(Shk) and U-band phot ometry-are known to be subop- 
timal for M dwarfs (e.g. Kiirste r et al.||2003 Dfaz et al. 
2007[ Gomes da Silva et clT 2011). Here, we have reexam- 


ined activity on Kapteyn’s star jsing the Ha and Na I D 
lines. We find that these lines offer a more significant de¬ 
tection of the P ~ 140-day periodicity observed by A14 
in Shk and FWHM. Based on the improved significance 
of the detection, and the mass and likely age for the star, 
we conclude the 143-day periodicity must be the stellar 
rotation period. This signal is strongly correlated with 
the proposed RV signal at P = 48 days. We conclude 
that the stellar rotation induces an RV signature which, 
when observed at the sparse cadence of the present RV 
set, creates the signal attributed to planet b. 
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2. DATA 

A14 claimed the detection of planets Kapteyn b and 
c based on RVs derived from 135 spectra obtained from 
the HARPS, HIRES, and PFS spectrographs. Our anal¬ 
ysis focused on the 95 HARPS spectra, which are all 
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Fig. 1. — Periodograms of the HARPS spectral activity tracers. The vertical scales are not identical. The peak near P = 140 days 
(highlighted) is present in every tracer except BIS. The dashed lines indicate the power required for a FAP of 1%. 


publicly availabl(|R We also considered the 20 (of 32) 
HIRES spectra ymere extracted spectra containing Ha 
were available in the Keck Observatory Archive. PFS 
spectra are not publicly available, and only account for 
8 of the available RVs. We adopted RVs from A14, 
which are derived using the HA RPS-TERRA algorithm 
(Anglada-Escude & Butler|2012 ) and a HIRES iodine-RV 
algorithm using these same spectra. 

We considered an array of stellar activity indicators 
for Kapteyn’s star. For the HARPS spectra, we mea- 
sure d the activity indices for the Na I D resonance lines 


(Jn: Gomes da Silva et al. 2011) and the Ha line (/h, 


Robertson et al.|2013 ), which have been shown to be par- 


ticularly sensitive to activity in M dwarf stars. We also 
consider the Shk values provided in A14. For HIRES 
spectra, only /hq is available. 

While the HARPS reduced data archive provides val¬ 
ues for the line-shape characteristic quantities Bisec¬ 
tor Inverse Slope (BIS) and full-width at half-maximum 
(FWHM), we have computed our own values for these ac¬ 
tivity tracers using the HARPS CCF (cross-correlation 
function) files, which are a product of the cross¬ 
correlation betwee n stellar spectra and a ta ilored numer¬ 
ical mask. As in Robertson et al. (2014), we summed 
the individual CCl's of orders in the wavelength range 
between ~5400 - 6200 A to obtain an “aver a.ge” CCF, 
ther i calcu lated BIS and FWHM according to [Wright et] 
al. (2013). Although the periodicites and correlations 
observed herein for BIS and FWHM have higher signifi- 

Data obtained from the ESO Science Archive Facility under 
request numbers 115242 and 149297 


cance when using our reduction, their qualitative behav¬ 
iors are identical for the original HARPS values. 

In addition to quantities derived from the HARPS 
spectra, we exa mined phot ometry of K apteyn’s star from 
the Hipparcos (|ESA||1997|) and ASAS (|Pojmanski||1997 ) 
archive s. Photometry of K apteyn’s star from KFL'I- 
South (Pepper et al. 2012) also exists, but the star’s 
proper motion is so high that it moves between pix¬ 
els throughout the observations, causing large offsets on 
the same timescale as the stellar rotation period (J. Ro¬ 
driguez, private comm.). We therefore do not incorpo¬ 
rate the KELT-S data into this analysis. 

We find that Shk and /d are clearly correlated with 
each other, and anticorrelated with Ha. There are three 
spectr£0for which the spectral indices do not follow these 
correlations. We excluded these spectra from our analy¬ 
sis, as they likely correspond either to flare events, other 
abnormal activity, or errors in the measurement of one 
or more lines. 

3. ANALYSIS 

3.1. Determining the stellar rotation period 

We began our analysis by examining the time series of 
the activity indices for any periodic signals which might 
leave an imprint on the RVs. In particular, for the peri¬ 
ods of the proposed planets, stellar rotation is the most 
likely source of astrophysical noise to create impostor 
signals. We computed the generalized Lom b-Scargle pe- 
riodogram (Zechmeister & Kiirster 2009) to search for 
periodicities in the activity tracers. 

5 BJD = 2454879.556, 2456424.446, 2456428.462 
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Fig. 2.— Top: Ilia measurements from HARPS (red) and HIRES (blue), with our fit to the 143-day stellar rotation signal (black). 
The inset (lower right) shows a detailed view of the observations conducted by A14. Lower left: /hq measurements folded to the rotation 
period. 


In Figure 0 we show that the Hct, Na I D, FWHM, 
and Shk inaices all show power near a period of 140 
days, as observed by A14 for FWHM and Shk- The sig¬ 
nal is strongest in Ina (Figure]^, where its periodogram 
power of 49.0 corresponds to a false-alarm probability 
FAP ~ 5 X 10“^® according to the approximation given 
in Equation 25 of [Zechmeister & Kiirster (2009). Pe- 
riodograms of the activity indices also show power at 
very long periods, potentially indicative of a long-term 
trend or curvature; indeed, visual inspection of the activ¬ 
ity series appears to show strong curvature as might be 
expected for a long-period cycle. However, as explained 
further below, we believe this is likely a spurious artifact 
caused by sparse sampling of the stellar rotation signal. 

We modeled sinusoids to the activity indices, using 140 
days as an initial guess for the period. For the combined 
HARPS/HIRES lua series, we found a best-fit period 
142.9 ±0.3 days with amplitude 0.0029 ± 0.0002 (= 4%). 
Our model to the data is shown in Eigurej^ Eits to other 
tracers produced very similar results. 

There are a number of outliers to this model, which 
may be caused by atmospheric distortion, errant data 
reduction, phase shifts in the rotation signal, or stochas¬ 
tic variability of the star. Of particular note are the 
two HARPS observations from October 2005 (phase ^ 
—0.22); the anomalous behavior is observed twice in 2 
nights, so the phenomenon is likely intrinsic to the star. 
The anticorrelation between Ha and the other tracers 
holds, precluding a flare. While these excursions are in¬ 
teresting, and may merit further study, the overall behav¬ 
ior of the stellar activity is well described by the 143-day 


rotation signal. 

It is necessary to justify our choice of a 143-day sinu¬ 
soid over a polynomial trend for our model to the activ¬ 
ity indices. First, we note that the periodogram power 
at long periods is driven almost entirely by the two ob¬ 
serving runs contributed by A14 in late 2013/early 2014, 
for which the activity indices are higher/lower (depend¬ 
ing on the tracer) than the earlier values. If-as sug¬ 
gested by our 143-day model (Figure inset)-those ob¬ 
servations coincidentally occurred near the same extreme 
phase of a large-amplitude rotation signal, it could create 
the appearance of a long-term trend. On a more physi- 
cally motivated level, w e note that at M* = 0.281Mq 
(Segransan et al. 2003), K apteyn’s star is below the 
threshold (iW* = 0.35Mq; |Chabrier & Baraffe||1997) 
where stars become fully c onvective. Long-term su rveys 


of activity in M stars (e.g. Robertson et al. 2013|) have 
confirmed that fully convective stars, which therefore 
lack the aSl magnetic dynamo believed to drive Solar- 
type activity cycles, tend not to exhibit long-period mag¬ 
netic cycles. 

Additionally, while the activity periodograms show 
many statistically significant peaks, we find that after 
fitting a 143-day sinusoid, the periodogram of the resid¬ 
uals shows no significant peaks, suggesting the power at 
other periods is simply aliasing from the rotation. When 
modeling other periods, especially long-term trends, con¬ 
siderable power remains in the residuals. We therefore 
conclude that a single 143-day sinusoid driven by the 
stellar rotation is the simplest, most physically plausible 
model for the observed variability. 
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Fig. 3. — RV as a function of activity for the HARPS-TERRA RVs and the predicted RVs from planets b and c. We do not show 
uncertainties for the model-predicted RVs. Our fits to correlated quantities are shown as red lines. We also show FWHM as a function of 
/hck- The strong anticorrelation confirms that stellar activity (traced by /hq:) alters the stellar line profiles (traced by FWHM), thereby 
shifting RVs. 


3.2. Impact of activity on RV 

A rotation period of 143 days for Kapteyn’s star is a 
source of concern for both of the purported planets in the 
system. The rotation period is uncomfortably close to 
the 121.5-day period of planet c, and almost exactly three 
times the period of planet b. We therefore considered 
possible correlations between the RVs and stellar activity. 

In Figure we show the HARPS RVs (which are the 
dominant source of information for the planet signals) 
plotted against /ho- Performing a linear regression, we 
find a Pearson correlation coefficient r = —0.38. The 
r-value suggests a probability P{r) = 9 x 10“^ that we 
would observe this level (or more) of correlation for 92 
measurements of these quantities if /ua and RV were in 
fact uncorrelated. 


Anglada-Escude & Tuomi (2015) have questioned our 
use of residual RVs when considering activity-RV corre¬ 
lations. While we disagree with their assertion that our 
statistical met hods invalidate the resu lts of our analysis 
of GJ 581 (see iRobertson et al.||2015|), we seek to avoid 
any statistical ambiguity here. Instead of modeling and 
subtracting planets or activity-RV correlations one at a 
time from the velocities and examining the residual ve¬ 
locities, we instead consider correlations between activ¬ 
ity and the modeled RV signatures of each planet as re¬ 
ported by A14. If the predicted behavior of either planet 
matches the observed activity variations, it must be dis¬ 
carded as an activity signal, regardless of the statistical 
methods used to identify it. 


Figure]^ shows the predicted RV contribution of each 
planet at the time of each observation as a function of 
the corresponding /hq measurement. For planet c, its 
RVs are uncorrelated with activity. On the other hand, 
the RVs of “planet b” are very strongly anticorrelated 
with /Ha- Our fit to the RVs of this signal as a function 
of /Ha is qualitatively similar to what we obtained using 
the original velocities, but with less scatter, suggesting 
the activity-RV correlation is mostly contained within 
the model to the 48-day signal, and the 121-day signal 
distorts the relationship. Fitting out stellar activity with 
a planetary orbit model at 48 days would explain the 
assertion in A14 that Kapteyn’s star is an especially low- 
jitter target. 

Other activity tracers show a qualitatively similar cor¬ 
relation with the predicted RVs of the 48-day signal. The 
only difference is that /na and BIS evince a negative cor¬ 
relation, while the other indicators are positively corre¬ 
lated. The fact that we observe simultaneous changes in 
both BIS and FWHM ensures that the line profiles are 
varying asymmetrically, and alter RVs. It is especially 
reassuring to see (Figure that FWHM is correlated 
with RV and the spectral line indices, confirming that 
activity is disorting the stellar line profiles-and, by ex¬ 
tension, the RVs-rather than some external phenomenon 
which affects both quantities similarly. In Table we 
show the Pearson correlation coefficients of all spectral 
tracers with each planet signal as well as the original 
data. The strong, consistent correlations offer conclusive 
evidence that the 48-day RV signal is not a planet, but 
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a harmonic of the 143-day stellar rotation period. 


3.3. Unresponsive activity tracers 

We find it remarkable that the stellar rotation signal 
and its correlation with RV appear so consistently across 
four different activity tracers. However, it is equally im¬ 
portant to point out that we see neither periodic vari¬ 
ability nor activity-RV correlations in the photometry 
from Hipparcos and AS AS, and do not recover the peri¬ 
odicity in BIS. The absence of features in these data are 
consistent with the findings of A14. A power spectrum 
of the ASAS photometry shows the peaks near 1 year 
and 1100 days mentioned in A14, but we find that those 
peaks match peaks from the window function (the peri- 
odogram of the sampling), confirming they are caused by 
poor sampling. 

The lack of a rotation signal in the photo metry is not 
surprising, as previous studies of M dwarfs (Kiirster et 
al. 2003 Robertson et al.||2014) have revealed magnetic 
phenomena in Ha which alter RVs without a correspond¬ 
ing photometric signature. Bisectors are more difficult to 
measure in M type stars due to the prodigious presence 
of line blends, and the lower HARPS signal -to-noise com¬ 
pared to FGK stars (Ba§turk et al. 2011), which likely 
explains both the lack ot signihcant periodicity and the 
less significant BIS-RV correlations. 

4. DISCUSSION 

What caused the appearance of a planet-like RV sig¬ 
nal at 48 days? The sampling cadence of the HARPS 
time series for Kapteyn’s star is insidious, allowing power 
from the rotation signal to leak to shorter periods. The 
unfortunate coincidence that the two 2013/2014 observ¬ 
ing runs both sampled the rotation signal near the same 
phase worsened matters, as it reduced the signihcance of 
the 140-day signal in both RV and the activity tracers. 
The presence of another RV signal at 121 days further 
obscures the velocity contribution from the stellar rota¬ 
tion. 


Bonffis et al. (2013) identified a 2-signal RV solution 
very similar to that of A14 using only the HARPS ob¬ 
servations from 2003-2009. Even for the complete data, 
most of the phase coverage of the 48-day signal comes 
from these observations; the new data from A14 only 
cover approximately a third of the 48-day phase curve. 
While the 2003-2009 data alone do not offer statistically 
signihcant detections of any periodicities, they are suffi¬ 
cient to see that the rotation signal and its harmonic at 
Troi/3-i.e. the period of “planet b”-are present during 
this epoch. In Figurej^ we show the periodogram of /hq 
from 2003-2009. The Prot/3 harmonic is the strongest 
peak in the power spectrum. In other words, for the 
hrst 30 HARPS observations, which are crucial to the 
appearance of the 48-day signal, stellar activity matches 
the periodicity of the “planet.” The deleterious effect of 
sampling in this data set illustrates an important require¬ 
ment for ultra-precise RV surveys searching for low-mass 
planets; a consistent, regular cadence is required to re¬ 
solve, model, and correct semi-periodic stellar signals in 
RV data. 

As to the origin of the 121-day period, we see little 
in our analysis to contradict the interpretation of that 
signal as a planet. Planets are known to exist around 
low-mass halo stars, as evidenced by the multi-planet 


^ , . 1 ^ .I 
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Fig. 4.— Periodogram of /jja HARPS observations from 
2003-2009. Peaks corresponding to the stellar rotation period {Prot 
= 143 days) and its harmonic at Prot/3—the period of “planet b” 
(Ply = 48 days)-are clearly distinguished. 


system around Kepler-444 (|Campante et al.||2015|), so it 
is reasonable that planets would exist around Kapteyn’s 
star. However, its period is worryingly close to the stel¬ 
lar rotation period and its aliases, and is equally vul¬ 
nerable to the sampling concerns outlined for the 48-day 
signal. Kapteyn c should probably be considered a planet 
candidate pending additional observations to confirm the 
physical origin of the 121-day signal. 

We direct particular attention to the stark contrast 
between the strengths of the activity-RV correlations in 
Iho! and In to those in the other spectral trace rs. This 
work, along s ide numerous previous studie s (e.g. Kiirster] 
et al. 1 120031 [Gomes da Silva et al.]|2011[ [Robertson et| 


^ 2013 2014) have shown that activity tracers redder 
than bAif-particularly Ha-often reveal activity signals 
in old M stars that do not appear in Ga H&K. Too fre¬ 
quently, these unrecognized activity signals have been in¬ 
terpreted instead as low-mass planets in the HZ. Enough 
examples have now been provided to conclude that an 
analysis of Shk alone is insufficient for ruling out false¬ 
positive planets around M dwarfs. It is our opinion that 
either highly confident knowledge of the stellar rotation 
period, or an analysis of /jiaj Id, and/or an equivalently 
suitable M dwarf activity tracer should be required for 
the confirmation of low-mass planets around M stars. 

Our measured 143-day rotation period for Kapteyn’s 
star is consistent with its previous descriptions as an old 
halo star. According to the age-rotation-Lx curve from 


Engle & Guinan (|2011|), this period suggests an age of 
approximately ItJ Gyr. We find the amplitude of the stel¬ 
lar rotation signal in the spectral tracers surprising for 
such an old, quiet star; with an amplitude of 4% in /hq,, 
it is more than twice the am plitude (1.7%) of the ro tation 
signal observed for GJ 581 ( Robertson et al.|[2014 ). Ad¬ 
ditional observations might determine whether the am¬ 
plitude varies, or is a perpetually strong signal. Also 
interesting is the fact that the extant data are consistent 
with the rotation signal staying in phase over the ~ 10- 
year span of the observations-more than 25 rotations- 
although again, the sampling makes it impossible to be 
sure this is the case. 


Kapteyn’s star joins Barnard’s star (Kiirster et al. 
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Tracer 

Variable 

Fit 

r 

P(r) 


Data 

Vr = 23±5 — 290±70 X Ilia 

-0.38 

9 X 10-® 

Ha 

Planet b (model) 

Vr = 42±3 — 540±40 X Ilia 

-0.79 

< 10-® 


Planet c (model) 

Vr = 7±4 — 90±60 X Ilia 

-0.16 

0.06 


Data 

Vr = —5±2 + 130±40 -^D 

0.37 

1 X lO-'^ 

Na I D 

Planet b (model) 

Vr = —12±1 -1- 260±20 X 7d 

0.80 

< 10"® 


Planet c (model) 

Vr = 0±1 + 10±30 X lu 

0.05 

0.32 


Data 

Vr = —110^40 + 30±io X FWHM 

0.31 

0.001 

FWHM (km s-i) 

Planet b (model) 

Vr = —250±20 + 77±7 X FWHM 

0.77 

< 10-® 


Planet c (model) 

Vr = —304^30 + 9±9 X FWHM 

0.11 

0.15 


Data 

Vr = —4±2 -1- 18±6 X Shk 

0.32 

9 X lO-'^ 

Ca H H&K 

Planet b (model) 

Vr = —6±i -1- 27±5 X Shk 

0.52 

5 X 10“® 


Planet c (model) 

Vr = 0±i -1- 2±5 X Shk 

0.04 

0.35 


Data 

Vr = 0±1 — 0.05±o.04 X BIS 

-0.16 

0.06 

BIS (m s-i) 

Planet b (model) 

Vr = 0.6±0.5 — 0.10±0.03 X BIS 

-0.31 

0.001 


Planet c (model) 

Vr = 0.3±o. 4 ~ 0.0±o.9 X BIS 

-0.10 

0.17 


TABLE 1 

Linear fits, Pearson correlation coefficients (r), and P-values for HARPS RVs (in m s“^) as a function of the spectral 

ACTIVITY tracers. HeRE, “DATA” REFERS TO THE ORIGINAL RVS, WHILE THE PLANET MODELS ARE THE PREDICTED RVS AT EACH 

observation for each PLANET ACCORDING TO A14. 


20031 and GJ 581 (Robertson et al. 20141 on a growing 
list of old M stars for which one or more activity signals 
are observed in Ha (and potentially other spectral trac¬ 
ers) and create RV shi fts, but have not bee n detected 
in existing photometry. Kiirster et al. (20031 attributed 
this behavior to localized regions on the stellar surface 
where strong magnetic fields inhibit convection, but do 
not create dark spots. This phenomenon deserves ad¬ 
ditional attention, as it now appears to be common in 
M dwarfs, and clearly poses a significant obstacle to RV 
detection of exoplanets in those systems. 


5. CONCLUSIONS 

We have analyzed the stellar activity of Kapteyn’s star 
using a variety of activity indicators. We found the star 
has a rotation period of approximately 143 days, and cre¬ 
ates a periodic signal in the Ina, Id, FWHM, and Shk 
indices. This signal is present in the RVs, and due to the 
sampling of the extant time series, creates a planet-like 
signal at 48 days (= Pr-ot/5)- The results of this anal¬ 
ysis, along with recent work on GJ 581 and GJ 667G 


highlight a number of precautions which must be taken 
to avoid false positives when searching for low-mass M 
dwarf planets with RV. For Kapteyn’s star, the signal of 
the claimed planet “c” does not seem to be correlated 
with stellar activity indicators, and we recommend con¬ 
tinued observations to confirm the planetary nature of 
this signal. 
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